Finite size effects in nonequilibrium wetting 



A C Barato 

The Abdus Salam International Centre for Theoretical Physics 
Trieste 34014, Italy 

E-mail: acardoso@ictp.it 



Abstract. Models with a nonequilibrium wetting transition display a transition also 
in finite systems. This is different from nonequilibrium phase transitions into an 
absorbing state, where the stationary state is the absorbing one for any value of the 
control parameter in a finite system. In this paper, we study what kind of transition 
takes place in finite systems of nonequilibrium wetting models. By solving exactly 
a microscopic model with three and four sites and performing numerical simulations 
we show that the phase transition taking place in a finite system is characterized by 
the average interface height performing a random walk at criticality and does not 
discriminate between the bounded-KPZ classes and the bounded-EW class. We also 
study the finite size scaling of the bKPZ universality classes, showing that it presents 
peculiar features in comparison with other universality classes of nonequilibrium phase 
transitions. 



1. Introduction 

It is well-known that phase transitions can only occur in systems with infinitely many 
degrees of freedom. For example, the Ising model in a finite volume does not display a 
transition, however using finite size scaling theory [1] one can obtain properties of the 
system in the thermodynamic limit by studying finite systems. A similar situation is 
observed for models in the so-called genuine nonequilibrium universality classes [2], as 
is the case of direct percolation (DP) [3J. 

On the other hand, a d— dimensional (with d > 2) layer Ising model infinite in d — 1 
dimensions and finite with size L in the remaining dimension displays a transition also 
for finite L. However, the transition for finite L is different from the transition in the 
thermodynamic limit: In the first case we have the critical behavior of the Ising model in 
(d— l)-dimensions and in the second in <i-dimensions. Similarly, nonequilibrium wetting 
finite systems also have a phase transition. Here the wetting transition is defined as an 
binding-unbinding transition controlled by a growth rate. The reason is that both phases 
are basically determined by the sign of the velocity of a free interface, giving a well- 
defined (but size- dependent) transition point, even for finite systems. This circumstance 
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has been pointed out some time ago by Munoz [I]. Nevertheless, the properties of this 
transition and its implications in finite size scaling were never studied in detail and this 
is the purpose of the present paper. 

Nonequilibrium wetting is a very rich and interesting example of nonequilibrium 
critical phenomena. It basically corresponds to the study of the Kardar-Parisi-Zhang 
(KPZ) equation [5] added with a soft-wall potential ( see [6HH] for reviews). The KPZ 
equation defines a robust universality class of nonequilibrium growing free interfaces 
[9j[10]. It is a Langevin equation for the interface height h(x,t), where x gives the 
position on a d- dimensional interface and t stands for time. The KPZ equation added 
with a term accounting for the presence of a soft-wall was introduced by Tu et al. [TT] 
and reads 



where C( x ?^) is a gaussian white noise and V(h) is a soft-wall potential given by 



The Laplacian term aV 2 /i(x, t) is related to surface tension, a is the average interface 
velocity at zero slope and the nonlinear term A(V/i(x, t)) 2 is the lowest order term that 
breaks the up-down symmetry P}[TU]. 

The wetting transition can be described as follows. If a is bigger than the critical 
value a c the interface will grow like a free KPZ interface, and the term exp(— h) becomes 
irrelevant after some transient. For a < a c the interface stays bounded to zero height 
and does not propagate. The order parameter for this transition is the ensemble average 
of exp(-h), which is zero for a > a c and non-zero for a < a c . 

For A = we have equilibrium wetting and equation ([TJ becomes the Edwards- 
Wilkinson (EW) p2| equation added with the soft-wall potential. This equation 
has been previously introduced by Lipowsky [13] in order to study dynamics in 
equilibrium wetting (see [T3HT7] for reviews on equilibrium wetting). For A ^ we have 
nonequilibrium wetting and it turns out that the set of critical exponents characterizing 
the transition depends on the sign of A. Therefore, we shall consider three universality 
classes: the bounded-EW class (bEW) and the bounded-KPZ classes (bKPZ+ and 
bKPZ— ), where the sign refers to A. 

By performing a Cole-Hopf transformation n = exp(-h) in equation (pQ) the 
following multiplicative noise Langevin equation is obtained, 



dh(x,t) 
dt 




d 



V(h) + ffV 2 /j(x, t) + A(V/i(x, t)f + C(x, t), 



(1) 



V(h) = exp(-h). 



(2) 



<9n(x, t) 



an(x, t) + n(x, t) 2 + uV 2 n(x, t) + n(x, t)((x, t), 



(3) 



dt 



where A = —a < 0. The above equation is called the multiplicative noise 1 (MN1) 
equation [6]. Performing the same transformation in equation (JT]) but with A = a > 
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the MN2 equation is obtained ©. Clearly MN1 and bKPZ- (MN2 and bKPZ+) are the 
same universality class. 

The most robust universality class of nonequilibrium phase transitions is the DP 
universality class [3]. The DP Langevin equation is very similar to the MN1 equation, 
the difference resides in the term multiplying the noise. In the MN1 equation it is 
proportional to n(x,t), while in the DP equation it is proportional to n(x,tY^ 2 . In 
spite of the similarity of both Langevin equations, the bKPZ- (or MN1) and the DP 
universality classes are essentially different [I]. Finite systems in the DP universality 
class have no phase transition, the stationary state is the absorbing one independently 
of the value of the control parameter. On the other hand, finite systems in the bKPZ 
universality classes still display a transition. This was clearly demonstrated by Munoz [I] 
with the solution of the corresponding one- variable Fokker-Planck equations (obtained 
by eliminating the Laplacian term) for the DP and the MN1 classes, which roughly 
corresponds to a zero- dimensional limit. 

In this paper we study the critical behavior of the phase transition taking place in a 
finite system and the finite size scaling of the bKPZ universality classes. We consider the 
microscopic model introduced in [18], which is a restricted solid on solid (RSOS) model 
with a hard- wall at zero height. We solve this model exactly for small systems with three 
and four sites, going beyond the zero-dimensional approximation from jl]. As expected, 
we demonstrate that the critical exponents of the phase transition in finite systems 
are different from the critical exponents characterizing nonequilibrium wetting in the 
thermodynamic limit. Surprisingly, we find that the phase transition in finite-systems 
does not discriminate between the bEW, bKPZ+ and bKPZ- universality classes. 

Moreover, we study the finite size scaling of bKPZ universality classes. We show 
that the survival of a transition also in a finite system leads to finite size scaling that 
is different from the finite size effects observed in the other genuine nonequilibrium 
universality classes [2]. In order to perform numerical simulations we use the single-step 
(SS) model with a moving wall [2"T||2"2"]. In this model the critical point is known exactly, 
leading to clean numerical results. 

The organization of the paper is as follows. In the next section we introduce the 
models and discuss their phenomenology. In Sec. Ill we solve the RSOS model exactly 
for the three-site and four-site cases. In Sec. IV we analyze the finite size scaling of the 
bKPZ universality classes with numerical simulations. We conclude in Sec. V. 

2. Models definition and their critical behavior 

First we define the microscopic realization of the bKPZ equation introduced by 
Hinrichsen et al. [18]. The model is defined on a one-dimensional discrete lattice of 
size L with periodic boundary conditions. To each site i of the lattice a random variable 
hi is attached. It can take any integer value and is interpreted as the interface height. 
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Figure 1. Possible transitions in a certain interface configuration of the RSOS model 
with a wall. The crosses indicate that the transition cannot be carried out. 



The allowed configurations are constrained trough the RSOS restriction, |/ij±i — hi\ < 1, 
which generates a surface tension. Usually, the initial condition is a flat interface at 
zero height, i.e. hi — for all i. The model evolves by random-sequential updates 
and the possible transitions are: evaporation hi — > hi — 1 and depositions hi — > hi + 1. 
Only transitions that satisfy the RSOS restriction are carried out. The deposition rate 
(transition probability per unit of time) is q, and the evaporation rates are p and 1, 
with p being the evaporation rate at a plateau (see Fig. [TJ. The hard- wall is introduced 
by forbidding evaporation events at zero height, in this way negative heights are not 
allowed. 

In order to explain the wetting transition within this model we first consider a 
free interface (no wall at zero height). In this case the interface will propagate with a 
constant velocity after some transient. For a fixed value of p, if q is bigger (smaller) than 
a certain value q c , the interface propagates in the positive (negative) height direction. 
For q = q c the interface velocity is zero. Now consider the wall is present. For q > q c 
it has no influence because the interface moves in the positive height direction away 
from the wall. On the other hand, for q < q c the situation is completely different: since 
the interface cannot propagate in the negative height direction, it stays bounded to the 
wall. Therefore, at q = q c that is a wetting transition, with q > q c corresponding to the 
the wet (or moving) phase and q < q c corresponding to the bound phase. In Fig. [2] we 
show the phase diagram of the model obtained from numerical simulations. We note 
that the critical line depends on the system size and we are denoting by q c the critical 
line in the limit L — > oo. When considering the critical point of a finite system we write 

It can be shown that the parameter A in the KPZ equation can be obtained within 
a microscopic model by calculating how the interface velocity vary with the interface 
slope [9]. For the RSOS model the slope is introduced by changing the RSOS condition 
between sites 1 and L in the following way. In order to have an interface slope m = M/L 
the RSOS condition between the boundary sites is changed to hi — hi = M, M+l, M—l. 
In Fig. [5] we also display the line A = obtained from numerical simulations, where 
above (below) this line A < (A > 0). Therefore, for < p < 1 the phase transition 
is in the bKPZ- universality class and for p > 1 in the bKPZ+ universality class. At 
p = 1 and q < 1 detailed balance holds (see [18]) and we have equilibrium wetting, i.e. 
the phase transition is in the bEW universality class. 

The order parameter of the phase transition is the density of sites at zero height 
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Figure 2. Phase diagram of the RSOS model. We show the lines v = and A = 
which corresponds to L — > oo and are obtained from numerical simulations. Moreover 
the figure shows the exact critical lines for L = 3 and L = 4, which correspond to 
equation (fT4]) and (fT9j) respectively. 



Po, defined by 

L 

Pb = <£<W>, (4) 

8=1 

where the brackets indicate an ensemble average. In the stationary state it is zero for 
q > q c and positive for q < q c . The critical exponent /3 is related to how the stationary 
value of the order parameter goes to zero as criticality is approached from below, i.e., 

Po ~ (5) 

where is the stationary value of the density of sites at zero height, A = q c — q, and 
the above relation is valid for A small and positive. At the critical point the order 
parameter as a function of time goes to zero following a power law with exponent 9, i.e. 

Po(t) ~ t- e . (6) 

The critical exponents defined by the divergence of the temporal correlation length Q\ 
and the spatial correlation length (± are 

C||~A^ii and C±~A _l/± , (7) 

which are valid for A small and positive. 

The KPZ universality class is related to the invariance of scale of the interface 
roughness pjITU] . It turns out that the critical exponents of the bKPZ universality classes 
are associated with the scaling exponents of the KPZ universality class [20] . Hence, it is 
known that, in one dimension, u\\ =3/2 and v^ = \ for the bKPZ universality classes and 
z/|| = 2/3 and u± = 4/3 for the bEW class [20]. From dimensional analysis one can see 
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that 9 = /3/vn, and, therefore, with the introduction of the wall only one new exponent 
arises (which can be 9 or /?). For the bEW this exponent can be determined exactly 
using the transfer matrix method and its value is (3 = 1 (or 9 = (3/v\\ = 3/4) [18J. 
For the bKPZ universality classes no exact solution is known and the exponents are 
determined with numerical simulations. To our knowledge the best numerical values for 
these exponents were obtained in [21] with time- dependent simulations of the SS model 
with a wall. The exponents are 9 = 1.184(10) (or (3 = 9v\\ = 1.776(15)) for the bKPZ- 
universality class and 9 = 0.228(5) (or (3 = 6v\\ = 0.342(8)) for the bKPZ+ universality 
class. 

Let us now define the notation for the distance from criticality and the critical 
exponents of a finite system. The distance from criticality will be denoted by A = 
q c (L) — q. An important result that we assume to be valid for the bKPZ universality 
classes and has been verified numerically in [20], for the RSOS model, is that 

\A-A\^L~ 1/v± . (8) 

The critical exponents for a finite system (3 and 9 are related to how the stationary 
value of po approaches zero as A — > and how po decays in time at the finite system 
critical point (A = 0), respectively. 

In this paper we also consider the SS model with a moving wall, which is more 
convenient to perform numerical simulations. In the following we describe this model 
briefly, for a full definition see [2"T| |2"2"] . Similarly to the RSOS case this model is defined 
on a one-dimensional discrete lattice with periodic boundary conditions, however in the 
SS model the height difference between two neighbors is \hi±\ — hi\ = 1. The possible 
transitions are: deposition in a local minima with probability s and evaporation in 
a local maxima with probability 1 — s. This model can be mapped onto the simple 
symmetric exclusion process and, therefore, its interface velocity is known and given by 

«(L) = (*-1/2)(1 + 1/L). (9) 

The moving wall is implemented by increasing (or decreasing) the height of the wall by 
one unit after a fixed time interval AT, which produces a wall velocity vw = 1/AT 
(vw = —1/AT). If the wall velocity is equal to v{L) we are at the critical point (of a 
finite system). We will consider the cases s = 1 and s = 0, which correspond to the 
bKPZ- and bKPZ+ universality classes, respectively [21]. All the critical exponents 
defined for the RSOS model are defined in the same way for the SS model, where the 
distances from criticality are given by A = \vw — v(oo)\ and A = \vw — v(L)\. Note that 
from we see that (JS|) is exact for the SS model with a moving wall. The advantages 
of performing numerical simulations with this model come from the fact that the critical 
point is known exactly. 

In the next section we solve the RSOS model exactly for the cases L = 3 and L = 4 
in order to make an exact analyses of the stationary state of a finite system. Before 
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ending this section we discuss the critical behavior of the RSOS model at p — 0. This is 
a special case where the phase transition is in the DP universality class [T2j. Comparing 
the transition at p = with the transition for p ^ provides a simple explanation for 
the essential differences between the DP universality class and the bKPZ universality 
classes. At p = once a layer is filled it is not possible to evaporate a particle on it 
anymore. Hence, in a finite system along p = for any q ^ the interface will have a 
positive velocity, since eventually a layer will be completed. Therefore, at p = a finite 
system has no phase transition. For p ^ particles can be evaporated from complete 
layers and because of that a phase transition is still present in a finite system. 



3. Exact solution for L = 3 and L = 4 

Let us consider the RSOS model with L = 3. We denote by P(h, h', h") the joint 
probability that site 1 has height h, site 2 has height h! and site 3 has height 
h" . For each height h there are seven possible configurations that satisfy the RSOS 
restriction. Because of translation invariance some of these configurations have the 
same probability, reducing the number of variables we have to deal with. For example, 
P(h, h, h+l) = P(h, h+1, h) = P(h+1, h, h). More precisely, we have three independent 
variables and they are: 

Xh = P(h, h, h) 
y h = P{h,h + l,h) 

z h = P(h+l,h,h+l). (10) 
Since there is a wall at zero height, Xh = Uh = z\i = for h < 0. 
The master equation for the three-site system reads 

= Uh + qzh-i ~[q + p(l - 6 h ,o)]xh 

j t Vh = ^ (qx h + 2z h - (2q + l)y h ^j 

^Zh = ~(2qy h +px h+1 - (q + 2)z h ^, (11) 

where the above equations are valid for h > 0. As an example, consider the first 
equation. The term comes from the transition (h, h + 1, h) — > (h, h, h), which is an 
evaporation taking place with rate 1, while the term qzh-i comes from the transition 
(h, h — 1, h) — > (h, h, h), which is a deposition taking place with rate q. On the other 
hand, the terms with minus sign come from the reversed transitions. The factor (1— 5h .o), 
also in the first equation, is related to the presence of the wall at zero height. The other 
two equations are derived in the same way. 

Because of the RSOS restriction, the one-site probability distribution and the 
two-site probability distributions Ph,h' are related by Ph = Ph,h + Ph,h+i + Ph,h-x- From 
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the relation between the two-site probability distributions and the three-site probability 
distributions we obtain 

Ph = %h + Vh-i + 2y h + 2^_i + z h . (12) 

The order parameter in the present notation is the probability of being at height zero 
Pq. The stationary state solution of the master equation ffTTl) gives 

p °" = w+w+w^ ~ ?3 ' (13) 

where the above solution is valid below the critical line, given by 

p = q 3 . (14) 

Therefore, for L = 3 the critical line is given by p = q 3 and /3 = 1 for any value of 
p. By integrating equations ( ITTj) numerically at the critical point we obtain 9 = 1/2, 
also independent of the value of p. Hence, the critical behavior for this finite system is 
different from the critical behavior obtained in the limit L — > oo. 

The case L = 3 is a very peculiar situation because one cannot introduce an interface 
slope and therefore the notion of A is lost. In this sense it is not surprising that the 
critical behavior is the same, independent of the value of p. For this reason we proceed 
solving the system for L — 4. 

In the case L = 4 there are six independent variables, and they are 
Xh = P{h, h, h, h) 
y h = P(h + l,h, h, h) 
z h = P(h + l,h+l,h+l,h) 
u h = P(h,h+l,h,h+l) 
v h = P(h,h + l,h+l t h) 

w h = P(h + l,h,h-l,h), (15) 
where Xh = yh = z h = Uh = v h = for h < and Wh = for h < 1. 
The master equation for L = 4 is given by 

^Xh = Vh + qzh-i ~[q + p(l - S hfi )]x h 

^Vh = q%h + u h + 2v h + qw h - [1 + p{\ - 5 KQ ) + 3q]y h 
d 

—z h = px h+1 + qu h + 2qv h + w h+1 - (p + 2 + 2q)z h 

^jU h = qy h + pz h - (1 + q)u h 
d 

—v h = 2qy h + 2z h - (2 + 2q)v h 
dt 

—w h = qz h _ 1 +py h -(q + l)w h . (16) 
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Figure 3. On the left panel we have the height profile obtained with numerical 
integration of the master equation (|16|) for L = 4 for the RSOS model at the critical 
point q = 2 and p = 6. On the right panel: the height profile for the SS model at the 
critical point in the bKPZ- case (s = 1) for L=32 obtained from numerical simulations 
(5 x 10 5 independent realizations). In both cases, the collapse of the data shows that 
the average height behaves like a random walk. 



The above equations are derived from the transition rules of the model, as we did for 
the three-site case. The one-site probability distribution as a function of the four-site 
probability distributions reads 

Ph = x h +3y h +y h -i+z h +3z h - 1 +u h +u h -i+2v h +2v h -i+2w h +w h -i+w h+ i.(17) 

The stationary state solution of equation f JT6|) gives 



p st = 2 + 3g(l + g)(2 + g 2 )+p(l + 3g + 3 g 2 ) , _ 4 

p(2 + p)[(2+p)(l + 4g + 6g 2 ) + 12g3( g+ lT P ^ P q 

which is valid below the critical line 



V 



-l + ^TTV- (19) 



Note that in both finite systems at p = we have q c = (see equations (fl~4l) and (fT9|) ). in 
agreement with the discussion in the previous section stating that at p = 0, for a finite 
system, there is no phase transition. As in the three-site case, the critical exponents 
are again independent of p, and given by $ — 1 and = 1/2 (obtained from numerical 
integration of (TT6j) at criticality). 

More precisely, by integrating the master equation numerically at criticality (for 
L = 3 and L = 4) we observe that the height profile is a gaussian scaling with t 1 / 2 : 
the average interface height performs a random walk with a reflecting boundary (see 
Fig. [3]). We also observed numerically (with Monte Carlo simulations) that the critical 
behavior that we obtained here with the exact solution of the cases L = 3 and L = 4 
holds for any finite system. As an example, in Fig. [3] we show the gaussian height profile 
for L = 32 obtained with numerical simulations of the SS model in the bKPZ- case at 
the critical point. Surprisingly, the transition in a finite system is the same for all the 
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three universality classes. We note that we can also write down the master equation 
for the case L = 5 (there are 10 independent variables), which is considerably more 
cumbersome. However, we could not find an analytical expression for the critical line in 
this case. 

We could consider a simpler case: a one site model, which would be a drifted random 
walk in the presence of a hard- wall (reflecting boundary). In this case, the drift is the 
control parameter and the critical point is at zero drift (where we have a reflecting 
boundary random walk). Therefore, already with this very simple approximation we 
have the critical behaviour of any finite system. Actually, this is the content of the 
one-component Langevin equation solved in [3], which is obtained by eliminating the 
Laplacian from the MN1 equation. In the height variable it corresponds to a random 
walk in the presence of a soft-wall potential. 

4. Finite size scaling of the bKPZ universality classes 

The purpose of this section is to study the finite size scaling of the bKPZ universality 
classes. Since there is a transition also in a finite system, the finite size scaling is 
different in comparison with the others genuine nonequilibrium universality classes [2J. 
For example, consider the DP universality class. In a finite system in the long time 
limit the system always goes to the absorbing state. The order parameter, in the active 
phase and near the critical point, as a function of time has initially a power law decay 
followed by a saturation at some value. If one considers times much bigger than the 
temporal correlation length (which depends on the system size) one sees a fast decay of 
the order parameter from the saturation value to zero. Clearly for the bKPZ universality 
classes the situation is completely different. We expect to see a crossover (in off-critical 
and time-dependent simulations) from the infinite system critical behavior to the finite 
system critical behavior. 

We show in Fig. H] time-dependent numerical simulations at the critical point (of 
the finite system) of the SS model for the bKPZ+ and the bKPZ- universality classes. 
From the data collapse in Fig. H] we see that the order parameter at criticality has the 
scaling form 



where f(x) is a scaling function that crosses over from f(x) ~ x~ e for small values of 
x to f(x) ~ x" 1 / 2 for bigger values of x. In other words, initially the order parameter 
decays with the exponent 9 that characterizes the critical behavior for L — > oo and 
for times considerably larger than the temporal correlation length the decay exponent 
crosses over to the random walk exponent 1/2. 

For the off-critical simulations the usual scaling form is given by 



p (t,L) = L~^f(tL- z ) 



(20) 




(21) 
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Figure 4. Scaled Order parameter po(t, V)L^' V± as function of the scaled time tL~ z 
at the finite system critical point obtained from numerical simulations of the SS model. 
The crossover from t~ 8 to t 1 / 2 is evident in the figure, for both universality classes. 




Figure 5. The scaled order parameter in the stationary state p^ t (A,L)L^^ u± as 
a function of the scaled distance from criticality AL 1 /^- 1 - obtained from numerical 
simulations of the SS model. 



where gt(x) is a scaling function. The scaling function is expected to follow gi(x) ~ x 13 
for small x. In Fig. E we plot p^(A,L)L^ x AL 1/u± for the bKPZ- and bKPZ+ 
universality classes obtained from numerical simulations of the SS model. We see a 
good data collapse and the scaling function having the expected form. 

Let us now analyse the order parameter in the stationary state as a function of the 
finite system distance from criticality. First we define the quantity 

ff(L,L) = rf{A,L). (22) 

From OH]) and fl2"T|) we expect the following scaling form 

pt(A,L) = L-^^ / (AL 1 /^), (23) 

where the scaling function <?/(x) should behave as gf(x) ~ x 13 (with j3 = 1) for small 
values of x. In Fig. [H]we plot this scaling function obtained from numerical simulations 
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Figure 6. The scaled order parameter in the stationary state, defined in (|22D . 
Po*(A, L)L I3 / V± - as a function of the scaled distance from the finite system critical 
point AL 1 /" 1 - obtained from numerical simulations of the SS model. 



of the SS model with a wall. For the bKPZ- case we see the scaling function approaching 
the expected linear form as we get closer to the critical point. For the bKPZ+ case we 
were not able to access regions close enough to criticality such that the form ~ x 
sets in, considering the last two points for each system size we get an exponent of the 
order 0.57, still far from 1. However, one can clearly see an increasing curvature as x 
gets smaller. 

We would like to stress three points. The first is that a similar situation for the 
off-critical finite size scaling has been considered in equilibrium in [23], that is: a system 
with an exponent /3 for L — > oo and a different exponent f3 for finite L, independent 
of the system size. The same scaling forms ( l2~Tj) and (|2"3"|) were demonstrated to hold 
in [23]. Secondly, we think that previous numerical results of the critical exponent /3 
for the bKPZ- class from off-critical simulations, as for example in [T5],|2~0"]. estimate a 
value smaller than 1.776(15) (obtained from time-dependent simulations in [21]) because 
the finite-size effects were not took into account properly. The third point is that we 
also performed simulations with the RSOS model and the results are in agreement with 
the results for the SS model, however in the second case we get much better numerical 
results. 



5. Conclusion 



As we cited in the introduction, it has been pointed out some time ago that the models in 
the bKPZ universality classes have a phase transition in finite systems [I], in opposition 
to the DP universality class. However, it was not clear what kind of transition takes 
place in a finite system. Here this issue was clarified: we showed that the phase transition 
in a finite system is characterized by a linear decay of the stationary value of the order 
parameter as criticality is approached and the average height of the interface performing 
a reflecting boundary random walk at the critical point. Moreover, we showed that the 
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transition in a finite system does not differentiate between the bKPZ+, bKPZ- and 
bEW universality classes. We have considered two microscopic realizations of the bKPZ 
equation: the RSOS and the SS model with a hard wall. However we expect this finite 
system phase transition to be a feature of the bKPZ universality classes. 

As another main result of this paper we have established the finite size scaling theory 
for the bKPZ universality classes, which is affected by the fact that finite systems also 
have a transition. For time-dependent simulations we showed that the scaling function 
(120]) has a crossover from the infinite system size decay exponent to the finite system 
size decay exponent. For the off-critical simulations we showed that the scaling function 
have different forms, depending on considering the distance from the infinite system 
critical point or the finite system critical point. In the first case the scaling function 
defined in (|2ip follows gi(x) ~ or for small x and in the second case the scaling function 
defined in (123|) follows gj{x) ~ x 13 for small x. 

Two remaining main open challenges in nonequilibrium wetting are the exact 
calculation of the new exponent, arising with the introduction of the wall, for the 
bKPZ classes and the experimental verification of the rich critical behavior that is 
predicted with the theory (see 0E] for a discussion). For the second challenge using 
the appropriate finite size scaling theory would be very important. 
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